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Abstract 

We have observed the primary electron spectrum from 30 GeV to 3 TeV 
using emulsion chambers flown by balloons at the top of the atmosphere, 
for the purpose of exploring the origin of cosmic rays in the Galaxy. The 
atmospheric gamma rays have been simultaneously observed in the 30 GeV 
~ 8 TeV energy range. In this paper, we estimate the atmospheric electron 
spectrum in the upper atmosphere (< 10 g/cm 2 ) from our observed gamma- 
ray spectrum using the electromagnetic shower theory in order to derive the 
primary cosmic-ray electron spectrum. The transport equations of the elec- 
tron and gamma-ray spectrum are analytically solved and the results are 
compared with those of Monte Carlo simulation (MC). Since we used the ob- 
served atmospheric gamma rays as the source of atmospheric electrons, our 
solutions are free from ambiguities on the primary cosmic-ray nuclear spectra 
and nuclear interaction models included in MC. In the energy range above 
several hundred GeV, the Dalitz electrons produced directly from neutral 
pions contribute around 10 percent to the whole atmospheric electron spec- 
trum at the depth of 4 g/cm 2 , which increases in importance at the higher 
altitude and cannot be ignored in TeV electron balloon observations. 
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nuclear interaction models 



Preprint submitted to Elsevier 



October 12, 2012 



1. Introduction 

The primary electrons above 30 GeV using balloon-borne emulsion cham- 
bers have been observed since 1960's [l|j2|j3[. The major purpose of primary 
electron observations is to investigate the transport and the origin of cosmic 
rays in the Galaxy, and so far one believes that the most likely sources of 
electrons are supernova remnants such as Vela SNR j^]. High-energy cosmic- 
ray electrons rapidly lose their energy because of the synchrotron radiation in 
the Galactic magnetic field and inverse- Compton scattering with the Galac- 
tic photons. These processes cause a steeper spectrum with the index of 
— (3.0 ~ 3.3). On the other hand, cosmic-ray protons have the spectral index 
of around —2.7. Therefore the atmospheric electrons increase the ratio to 
primary electrons because they have the same spectral index as the protons. 
In order to derive the primary electron spectrum in balloon experiments, it is 
required to subtract the atmospheric electron spectrum produced in residual 
atmosphere. 

Atmospheric electrons above 10 GeV are mostly due to the pair produc- 
tions of gamma rays produced in the hadronic interactions in the residual 
atmosphere. Those gamma rays are mainly due to the decays of 7T° mesons 
with smaller fractions from rj mesons, K mesons, etc. In this paper we 
estimate the atmospheric electron spectrum from the observed gamma-ray 
spectrum and show the contribution of the Dalitz electrons which are di- 
rectly produced electrons from ir° and r\ mesons, that is not through gamma 
rays. In addition, we consider the direct production of electrons from K 
and tt^ mesons, which contributes in the lower energy range below 100 GeV 
and increases with decreasing energy. The first precise estimates of atmo- 
spheric electron flux was calculated by Orth & Buffington j4|. They have 
given the atmospheric electrons in the 1-100 GeV energy range contributed 
from the decay of mesons produced by nuclear interactions. Their calcula- 
tions are based on MC techniques and they presented just a simple analytic 
model which summarize their MC results. On the other hand, we estimate 
the flux of atmospheric electrons from the observed gamma rays by using 
the electromagnetic shower theory. The results are compared to the MC re- 
sults assuming the primary hadron flux and nuclear interaction model, and 
discuss the behavior of both electron and gamma-ray spectra in the upper 
atmosphere. 
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2. Calculations 



We consider the transport equations of the electron and gamma-ray spec- 
trum above 10 GeV in the atmosphere (< 10g/cm 2 ) based on the electro- 
magnetic shower theory of Approximation A |5| which treats the behavior of 
electromagnetic interaction in the upper atmosphere. The atmospheric elec- 
trons are mostly produced by gamma rays, originated from n° decay mode, 
rj and K° decay mode, etc. The equation should include the production rate 
of atmospheric gamma rays from those mesons. We estimate the production 
rate from the atmospheric gamma-ray spectrum J^{E) observed simultane- 
ously in our primary electron experiments p: 



J 7 (£)(m 2 -s-sr -GeV)' 1 

= (1.12 ±0.13) x 10- 4 (£/100GeV)- 2 - 73±0 - 06 , (1) 

in the range from 30 GeV to 8 TeV, which is the vertical gamma-ray spectrum 
normalized at 4.0 g/cm 2 originating from hadronic interactions of primary 
protons and heavy nuclei with the residual atmosphere. The contribution 
from primary electrons is subtracted. 

2.1. Formulation and Solution 

At the depth t in the atmosphere, the number of electrons, 7r(E,t)dE, 
and gamma-rays, ^(E, t)dE, with energies between E and E + dE satisfy the 
following simultaneous equations. 

<ML ' t] = -A'ir(E,t)+B'rtE,t)+* D (E,t) (2) 



dt 
di(E,t) 



C'n(E,t)-a ol (E,t)+ 7ex (E,t) (3) 



dt 

As shown in the literature the operator A' represents the effect of the 
radiation loss process, B 1 the pair production, C the gamma ray contributed 
from electrons, and o"o the absorption by the pair production. The functional 



forms of these operators are shown in |Appendix A 



In eq. (J3J), we add the gamma-ray production rate from nuclear interac- 
tions as: 

7ex (E, t)dE = k g E~P- 1 eM~t/L)dE , (4) 

in which the spectral index = 1.73 is determined from eq. (JT|). The attenua- 
tion path length is assumed as L = 100 g/cm 2 [6], giving the uncertainty of a 
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few percent to the production rate at the depth of 10 g/cm 2 . The coefficient 
kg is estimated from the observed gamma ray spectrum of eq. (pQ). Gamma 
rays are mostly produced through the 7r° decay mode and other small con- 
tributions exist from the rj and K° decay mode, which give the fraction of 
0.16 and 0.03, respectively, to the n° decay mode for the spectral index of 
around 1.7 0]. 

In the eq. ([2]) we include the small contribution of electrons from the 
nuclear interaction of Dalitz decay mode, which is the second major mode 
of 7T° decay, producing one real gamma ray and one electron-positron pair, 
occupying 1.2% of ir° decay process 0. The production rate of Dalitz electron 
pairs is estimated as a fraction of the gamma ray production rate as follows: 

7r D (E,t) = a-j ex (E,t) . (5) 

The coefficient a is estimated from the decay rate of Dalitz electrons from the 
7T° and rj mesons, and the result is given by a = 4.9 x 10 -3 (see Appendix B). 



The initial condition of electron spectrum at the top of the atmosphere 
is the primary electron spectrum [l[ and that of the gamma-ray spectrum is 
the Galactic diffuse gamma rays, which is assumed to be zero. We have 

tt(£,0) = 1.6xl0- 4 ( %— )- 3 ' 3 (6) 

v ' ; v 100 GeV ; V ; 

(m 2 • sr • sec • GeV) -1 , 
= 0. 

The eq. (jSJ) gives the electron spectral index of a = 2.3. 

The transport equations are solved using the above conditions. The so- 
lution of electron spectrum at the depth t becomes 

n(E,t) = Tr(E,0){(t)+ lex (E,t)£(t) , (7) 
[(A: + a )e^ - (A 2 + a )e^) 



Ai — A2 

B(P) e^+T^-l e ( A2+ r)*-l„ 



Ai-A 2 L Ax + l/L A 2 + l/L 
a 

1 Ax + l/L 
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(A 2 + ao)e^ +1 / L ) t -(a -l/L) 

A 2 + 1/L J ' 

The coefficient ((t) represents the energy loss rate of electrons by bremsstrahrung 
in the atmosphere. The term j ex (E,t)£,(t) represents the atmospheric elec- 
tron spectrum. The first term of £(t) represents the pair production rate from 
gamma rays and the second term including a is the contribution from Dalitz 
decay mode. The atmospheric gamma-ray spectrum at the depth t g/cm 2 is 
given by 

7 (£,t) = ^,0)^(1) +^ ex (E,t) V2 (t) , (8) 

A\ — A2 

V2(t) = ^-^-[/(A 1 )-/(A 2 )] , 

fw = Y ^ T j Z {^+m)e {x+l/L)t 
-{m-\) + *-cme {x+llL)t -i)} , 

where r]\ represents the gamma-ray production rate from primary electrons 
and 7/2 represents the attenuation rate of gamma rays by pair production. 
The coefficient Ai(s), A2(s), A(s), B(s), C(s) and <r with s — a, (3 used 



above equations are well known formulas and given in Appendix A 



The approximate expression of atmospheric electron spectrum becomes 
f (t) ■ lex (E, t) ~ (^±t 2 + a ■ t) lex {E, t) , 



which corresponds to the estimate of atmospheric electron spectrum in the 
Orth&Buffington's paper jij. 

The production rate of gamma rays eq. (J4]) estimated from the eq. (151) 
and the eq. (p]) is given by 

% x (E,t) = 1.11 x l(r 3 (£/100 Gey)" 2 ' 73 exp ' 



100 g/cm 2 
m 2 • s • sr • GeV • r.l.) -1 , 



in which the radiation length of air is taken as 36.6 g/cm 2 Q. 
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2.2. Solution with Zenith Angle Distribution 

As cosmic rays come in all directions, the intensities in the previous sec- 
tion should be integrated with the acceptance angle of each detector. At 
the depth t , the electron differential spectrum j ob (E,t) and the gamma-ray 
spectrum r y ob (E,t) are derived as follows. 

j ob (E,e,t) (m 2 • sec • GeV- 1 ) 

f 9 t 

= 2ir tt(E, -)cos0sin0d0 

Jo cos 9 

= ir(E, 0)CK 0, t) + lex (E, t)t((3, 9, t) (9) 

lob (E,9,t) (m 2 • sec • GeV" 1 ) 
t 



= 2ir [ 7 (£, 
Jo 



cos 9 sin 9d9 



cos 9' 

= tt(E, 0)m(a, 9, t) + lex (E, t)r] 2 {p, 9, t) 

The coefficients of each term are calculated from the following expressions, 

((a,9,t) = O 1 (a,9,t) + a O 2 (a,9,t) 
t((3,9,t) = B(0)O 3 (P,O,t) 

+a-{O 2 (P,9,t) + a O 3 (P,9,t)} 
77i(a,0,t) = C(a)0 2 (a,9,t) 
V2(P,6,t) = 2 (f3,9,t) 

+(A(/3) + a-C(f3))0 3 (f3 : 9,t) 

and the series of Oi(s, 0, t) with I = 1, 2, 3 are given by 

oo 

O^s, 0, t) = 2vr K(s)T n „ 1+l (9, t) 

n=0 

1 _ ms 2 f) 

T = - , T 1= t(l-cos0), 

t 2 1 f 3 1 

T 2 = ^og(— T s = l -(—-l), 
2! cos0 3! cos0 

t n 1 1 
... t = - —( - 1) ■■■ 

n n\n-2 K cos n - 2 9 ' 

\ \ \ n + 1 \ n + 1 1 

^0 = T 7~, • • • , A n = h A n _!—— , • • • 

Ai — A 2 A\ — A 2 L 
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Figure 1: Atmospheric gamma-ray spectra J g (E) of eq. (jHJ are compared with the MC 
results. The solution at 4 g/cm 2 is equal to eq. multiplied by E 3 . The MC data divided 
by 0.8 are also shown and good agreement with the analytic solution. 



3. Results and Discussions 

The analytic solutions obtained here are compared with the results of 
Monte Carlo simulation(MC). We use the COSMOS simulation code js|] in 
which the Dpmiet3 jsj is adopted as the hadronic interaction model and the 
BESS-TeV [10] 11] and RUN JOB 12| data for the input proton, helium and 
CNO spectrum at the top of atmosphere. We have confirmed the consis- 
tency of both methods for the cross section using the fundamental processes 
as shown in Appendix A In addition, we consider atmospheric electrons 



directly originating from K and tt^ mesons, which mainly contribute be- 



low lOOGeV. The decay modes of n — > 



e ± vv, 



Kl 



->• 



7T ± e =F ZA 



K- 
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— > 7r e v and so on, are calculated using the recent accelerator data 
and those contributions are treated as the ratio to electrons produced 



by gamma rays. 

3.1. Atmospheric Electron and Gamma-ray Spectrum 

The analytic solution of atmospheric gamma-ray spectrum given by 



J g (E,t) = lex (E,t)r] 2 (t) 



(10) 
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Figure 2: Atmospheric electron spectra are shown at the depth 4, 10 g/cm 2 . The analytic 
estimation from all mesons (eq. (ITT1) ') and the analytic solution £(t)~f ex (E,t) estimated 
from gamma-ray spectrum with and without Dalitz decay mode are shown with MC results 
divided by 0.80. 

in eq. (|HJ) and is shown in Figure [I] with the MC results at the depth of 
4, 10 g/cm 2 . Figure [1] shows that the gamma-ray spectrum of MC is 20% 
lower than the analytic solution, which has been indicated in our previous 
paper 0: the proton spectrum deconvolved from our observed gamma-ray 
spectrum assuming the nuclear interaction model of Dpmjet3 is 20% higher 
than observed primary proton data. Thus the discrepancy of gamma-ray 
flux in this paper is consistent with the difference between deconvolved and 
observed proton flux as described in the previous paper 

The analytic estimation of atmospheric electron spectrum at the depth t 
is given by 

j sec (E,t) = (l + C)^ ex (E,t) (11) 

as shown in eq. ([7]). The fraction C represents the contribution from direct 
productions from K and 7r =l= mesons and is shown as the increase of the thin 
line(all contributions) to the thick line(from gamma rays) in Figure [3j In 
the fraction C, the K\ mode are mostly dominant as indicated by Orth & 
Buffington [ijj. The contribution from K mesons estimated from the recent 
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Figure 3: The ratio of atmospheric electrons to gamma-rays. The analytic estimation 
from all mesons and the analytic solution from gamma-ray spectrum with and without 
Dalitz mode are compared with Monte Carlo simulation. The ratio r represents £(t)/rj2(t) 
and the fraction C represents the contribution from direct productions from K and tt^ 



accelerator data of K/n ratio 13] is a factor of two smaller than their result. 

The atmospheric electron intensity is approximately proportional to t 2 . 
The comparison with MC is shown in Figure [2} in which the spectrum of 
MC is divided by 0.80 to correct the discrepancy of gamma-ray spectrum. 
In the energy range above several hundred GeV, we notice that the Dalitz 
decay mode contributes more than 10 % to the ratio at high altitude above 
4 g/cm 2 , and the estimated spectrum is good agreement with the corrected 
MC simulation in all energy range. The atmospheric electrons become a 
larger fraction of the observed electrons in the higher energy region because 
of the harder spectrum than the primary electron one. 

We investigate the consistency between the analytic estimation and MC 
by comparing the ratio of atmospheric electrons to gamma rays, which is able 
to exclude the discrepancy of gamma ray spectrum. The ratios compared 
at the depth t = 4, 10 g/cm 2 are shown in Figure [3j which shows good 
agreement between the analytic estimation and MC simulation. The Dalitz 
contribution has constant value and does not vary with the depth t, since 
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Figure 4: The correction term of atmospheric electrons of eq. (fT5|) with acceptance angle 
of 30°, 45°, 60° and vertical for detectors at 4, 10 g/cm 2 , and compared with the result 
of primary electron spectrum. 



Dalitz electrons approximately increase in proportional to gamma rays and 
t, while electrons produced by gamma rays are proportional to t 2 . It means 
that the Dalitz mode plays a more important role at the higher altitude 
above several hundreds GeV. On the other hand, below lOOGeV the direct 
productions from K and n^ 1 mesons become larger than that from the Dalitz 
mode. 

3.2. Observation and Correction 

The primary electron spectrum is obtained from 

tt(^,0)= - lex {E,t) — m -sr- sec -GeV , 12 

C{a,e,t) C(a,0,t) 

using the eq. fl9]), in which l/£ is the correction term of the energy reduction 
by bremsstrahrung. The second term is the correction term of atmospheric 
electrons which are observed at the depth t within zenith angle 9 and cor- 
rected by energy reduction to the top of the atmosphere. The term £/£ is 
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Figure 5: The primary plus atmospheric electron spectrum (thin line) at each depth 
above 10 g/cm 2 including the extrapolation above 3 TeV. The primary electron spectrum 
(thick line) with the simple power law is estimated form the recent result below 3 TeV 0]. 
The atmospheric electron spectra (dashed line) are given by eq. (|13|) . 
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approximately proportional to t 2 . The eq. ( 1T2"j) indicates that all the observed 
electrons are firstly corrected by bremsstrahrung loss energy, and then, the 
atmospheric electrons are statistically subtracted. 

Finally the atmospheric electron spectrum corrected at the top of the 
atmosphere is given by 

J sec (E,t) = (1 + G) lex {E,t) 7^44 K • sr • sec • GeV ) _1 • ( 13 ) 

Figure H] indicates the correction term of atmospheric electron spectra of 
eq. ( JT3l) for the different acceptance angle of each detector at 4, 10 g/cm 2 , 
compared with the primary electron spectral data. The fraction C practically 
has no influence on the correction of the primary electron spectrum. Further 
in Figure EJ we show the primary electron spectrum, the atmospheric electron 
spectra of eq. ( Tl3l) and total spectra at each depth above 10 g/cm 2 . It is 
shown that the balloon experiment for electrons above 1 TeV requires the 
observation at a high altitude, since correction terms dominate the data at 
10 g/cm 2 . 



4. Conclusions 

In the balloon observations of cosmic-ray electrons, the atmospheric elec- 
trons are inevitably included in the observed spectrum and necessary to be 
subtracted to obtain the primary electron spectrum. The atmospheric elec- 
trons are mostly produced by the atmospheric gamma-ray conversion, and 
above several hundred GeV, around 10% of them at the depth of 4 g/cm 2 
are produced by the Dalitz decay mode. We have treated these processes 
analytically and obtained accurate solutions. The result shows that the at- 
mospheric electron spectrum obtained by MC is 20% lower than that of 
the analytic estimation. The discrepancy comes from the difference of at- 
mospheric gamma-ray spectrum, since the ratio of atmospheric electrons to 
gamma rays shows a good agreement. The difference of the gamma-ray spec- 
trum comes from the uncertainties of the primary cosmic-ray hadron spectra 
and the nuclear interaction models used in the MC calculations, together with 
the statistical error of 10 % included in the observed gamma-ray spectrum 
of eq. ([1]). The Dalitz mode gives a constant contribution to the atmospheric 
electron/gamma-ray ratio and becomes particularly important at the higher 
altitude above 4 g/cm 2 . On the other hand, in the lower energy range be- 
low 100 GeV the direct contributions of electrons from K and mesons 
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become larger than the Dalitz contribution, although atmospheric electrons 
below 100 GeV have almost no influence on the correction of the primary 
electron spectum. 

The astrophysical significance to search for the cosmic ray sources by pri- 
mary electrons become more important at higher energy regions, however, 
the contribution of the background atmospheric electrons increases seriously 
in the balloon observations in TeV region. Thus the observation of TeV 
electrons requires the higher altitude, and more precise estimates of the at- 
mospheric electrons are required than that of several hundred GeV region. 
Although these conditions are getting more difficult at higher energy region, 
the experiments in TeV region are the most important to search for cosmic- 
ray sources j^j. 
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Appendix A. Fundamental Parameters 

We calculate the fundamental process and compare the result with that 
of MC to confirm the consistency of cross section. When gamma rays or 
electrons pass through the thickness t of the air, the number of electrons 
tt(E, t)dE and gamma rays j(E, t)dE with energies between E and E + dE 
satisfy the transport equations as: 

ME ' t] = -A'*{E,t) + B' 1 {E,t) 



dt 
dj(E,t) 

dt 



C'n(E,t)-a ol (E,t) 



The physical meaning of each parameter A', B\ C and is explained in 
the section [27X1 in this text. The functional forms of the operators A', B', C, 
o"o arid A1.2 used in the equations are given in several references (e.g. [5j) as 

A(s) = 1.3603-^ lnlYs + 2) - £ 7 - 0.07513 , 

v ; ds K J (s + l)(s + 2) 



B(s) 



1.3603 



'8 + 1 (s + 2)(s + 3)' 



1 1.3603 

C ( s ) = + 



s + 2 s(s + l) ' 
a = 0.7733 , 

Ai, a (a) = \l-(Ms) + *o) ± {(A(s) + a ) 2 

-A(A(s)a - B(s)C(s))}^} , 

where s + 1 corresponds to a power-law index of the solutions. The numerical 
values of A(s), B(s) and C(s) used in this paper are listed as 

A(2.30) = 1.6743 , 5(2.30) = 0.4867 , (7(2.30) = 0.4118 , 

A(1.73) = 1.4270 , 5(1.73) = 0.5784 , (7(1.73) = 0.5561 . 

We calculate the gamma ray conversion process and confirm the agree- 
ment between analytical method and MC. If the initial condition is given 
by 

tt(E,0) = 0, 
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Table A.l: The number of Electrons produced by gamma rays at 100 GeV, which is 
divided by the initial number of k g . 



Depth (t) 


0.4 g/cm 2 


4.0 g/cm 2 


Analytic Solution 


6.229xl(T 3 


5.596xl0~ 2 


MC 


(6.25 ± 0.04) x TO" 3 


(5.61 ± 0.04)xl0" 2 



with (3 = 1.73, the solution of produced electron and gamma ray spectrum 
at the depth t becomes 

n(E,t) = ^^5(/3)[e^-en 

Ai — A2 
Al — A2 

+A ie Al * - A 2 e A2t ] . 

In Table \KA\ the number of electrons n(E,t) produced by gamma ray con- 
version are shown at the depths of 0.4 g/cm 2 and 4.0 g/cm 2 and compared 
with the results of MC. These values directly represent the cross section of 
pair production and one can find that the cross section agrees with each other 
within ±0.5 %. 



Appendix B. Dalitz electron contribution 

We include the contribution of Dalitz electrons represented by a fraction 
of gamma-ray production rate (eq fl5])). The coefficient a is estimated from 
the sum of Dalitz electrons of 7r°, r] and K° decays. 

In the 7T° decay, the branching ratio of the first mode is given by B m = 
B(n° 77) = (98.798 ± 0.032)% and that of the second mode is B D = 
B(n° ->- e+e~7) = (1.198 ± 0.033)% 0- We require the the ratio of Dalitz 
electrons to gamma rays from the 7r° meson spectrum oc E" 13 " 1 (j3 = 1.73), 
and it becomes B D /(2B m )/(2/(/3+l)) = 4.4x 10~ 3 , if the distribution of pair 
electrons is uniform. Precise estimate gives a somewhat different distribution 
of the pair production with two peaks, so that the value is slightly corrected 
and becomes 4.6 x 10 -3 . 

The ratio of Dalitz electrons to gamma rays for each decay is listed in 
Table IB.2I The major contribution of Dalitz electrons comes from 7r° decay. 
The rj decay gives the small contribution, and the K° decay can be neglected. 
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Table B.2: Dalitz contributions of tt°, r\ and K° decays 





TV® 


77 


K° 


Gamma-ray (77) mode 


98.8 % 


39.3 % 


98.8 % 


Dalitz (e + e~7) mode 


1.2 % 


0.68 % 


1.2 % 


Dalitz electron / 7-ray (r) 


4.6 x 1(T 3 


6.6 x 1(T 3 


4.6 x 10~ 3 


Fraction of 7-ray production (p) 


1/1.19 


0.16/1.19 


0.03/1.19 


Ratio in upper atmos. (= r x p) 


3.9 x 1(T 3 


0.9 x 10~ 3 


0.1 x 10~ 3 



The net fraction of Dalitz decay electrons to the observed gamma rays are 
given by 4.9 x 10~ 3 . 
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